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Abstract. Given a convex body C C containing the origin in its interior and 
a real number r > 1 we seek to construct a polytope P C C with as few vertices 
as possible such that C C tP. Our construction is nearly optimal for a wide range 
of d and r. In particular, we prove that if C = — C then for any 1 > e > and 
t = 1 + e one can choose P having roughly e~ d / 2 vertices and for r = \fed one can 
choose P having roughly d 1//e vertices. Similarly, we prove that if C C K d is a convex 
body such that — C C fJ.C for some /u > 1 then one can choose P having roughly 
(O + 1)/(t - l)) d/2 vertices provided (r - + !)<!. 



1. Introduction and main results 

We discuss how well convex bodies (compact convex sets with non-empty inte- 
rior) can be approximated by polytopes (convex hulls of finite sets of points). There 
is, of course, a vast literature on the topic, as there are many different notions of 
approximation, see surveys [G93a] and [Br07]. Our setup is as follows. Let C C M d 
be a convex body containing the origin in its interior. We seek to construct a 
polytope P cM. d with as few vertices as possible, so that 

P C C C rP 

for some given r > 1. 

Our first main result concerns symmetric convex bodies C for which C = —C 
and r measures the Banach-Mazur distance between P and C. 
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(1.1) Theorem. Let k and d be positive integers and let r > 1 be a real number 
such that 

(r _ v^TT)* + ( T + v^TT)' > 6 ( d + fc ) 1/2 . 

T/ien /or any symmetric convex body C C IR d t/iere is a symmetric polytope P C M d 
with at most 

'd + k? 
k 
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vertices such that 



P C C C tP. 



'd + P 

In fact, we can replace I ) throughout the statement of Theorem 1.1 by a 



slightly smaller number 



k 



a-") ^^E( d+ ::L" 2m )- 

m=0 V 7 

For example, taking d = 20 and /c = 3 we conclude that any 20-dimensional sym- 
metric convex body can be approximated within a factor of r = 3.18 by a symmetric 
polytope with at most 12, 480 vertices. 

Taking r in Theorem 1.1 arbitrarily close to 1, we obtain the following corollary. 

(1.2) Corollary. For any 

7 > w 0.48 

4v/2 

t/iere exists eo = £0(7) > such that for any < e < eo and for any symmetric 
convex body C C W 1 , there is a symmetric polytope Pel" 1 with at most 



Ve e 
vertices such that 

P C C C (l + e)P. 

The well-known volumetric argument (see, for example, Lemma 4.10 of [Pi89]) 
produces polytopes with roughly (3/e) d vertices which approximate a given sym- 
metric d-dimensional convex body within a factor of 1 + e. Hence for small e > the 
estimate of Corollary 1.2 gives us roughly the square root of the number of vertices 
required by the volumetric bound. It follows from results of Dudley [Du74] and 
also from results of Bronshtein and Ivanov [BI75] that in any dimension d one can 
construct a polytope P with not more than , y(d)e~( d ~ 1 ^ 2 vertices approximating 
a given symmetric convex body C C M. d within a factor of 1 + e, with 7(d) of the 
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order of d d l 2 . If the boundary of C is C 2 -smooth then for all sufficiently small 
< e < 6q(C) one can obtain an approximating symmetric polytope with at most 
(7 1 'ty ' '/ 2 vertices for some absolute constant 7 > 0, and the dependence on e 
cannot be made better [G93b], [BoOO] (note that the upper bound for e depends on 
the convex body C). The estimate of Corollary 1.2 is the first bound improving the 
volumetric bound uniformly over all symmetric convex bodies C of all dimensions 
d. 

Next, we consider approximations for which we want to keep the number of 
vertices of the polytope polynomial in the dimension of the ambient space. 

(1.3) Corollary. For any 

7 > y~ « 0.82 

there is a positive integer ko = ^0(7) such that for any k > ko and for any symmetric 
convex body C C M, d of a sufficiently large dimension d > do(k) there is a symmetric 
polytope P C M, d with at most 

'd + k'' 
k 

vertices such that 

P C C C 7 

A simple computation shows that if C is the <i-dimensional Euclidean ball and P 
has at most d k vertices for some fixed k, then P cannot approximate C better than 

r~d~ 

within a factor of r = 7W as d grows, where 7 > is an absolute constant. 

V Kind 

Finally, we consider approximations of not necessarily symmetric convex bodies. 
We prove the following main result, generalizing Theorem 1.1. The quality of 
approximation depends on the symmetry coefficient of the convex body C, that 
is on the smallest /j, > 1 such that — C C jiC (recall that the convex bodies we 
consider contain the origin in their interior). 

(1.4) Theorem. Let d and k be positive integers. For r, fi> 1 let us define 

A = A(r, /i) = — - — r + ^—Jl > 1. 

If 

(A-^T) fc + (A + ^T) fc > ^ + k k ) 1/2 

then for any convex body C C M d containing the origin in its interior and such that 
—C C fJ,C there is a polytope Pcl d with at most 
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d + k 

k 
3 



vertices such that 



P C C C tP. 



We also obtain the following extension of Corollary 1.2. 

(1.5) Corollary. 

(1) For t, fj, > 1 let us define 

S = 5{t, (jl) = 



fi+1 

For any 

7 > « 0.48 

there exists So = 5o(t) > swc/i t/iat as /ono as <5(r, //) < So, for any convex 
body C C M. d such that —C C jiC there exists a polytope with at most 



vertices such that 
(2) For any 



P C C C tP. 



7 > | w 0.34 



t/iere exists eo = £(7) > such that for any < e < eo and for any convex 
body C C W 1 such that — C C fxC for some fj, > 1 t/iere exists a polytope 
P GR d with at most 



vertices such that 

P C C C (l + e)P. 



As a function of the symmetry coefficient //, the number of vertices of P grows 
roughly as [i d ^ as long as the ratio t//j, is small enough. It follows from results of 
Gruber [G93b] that if the boundary of C is C 2 -smooth then for all sufficiently small 

< e < €q(C) one can construct a polytope P with not more than // d / 2 (7/e)( d-1 - ) / 2 
vertices for some absolute constant 7 which approximates C within a factor of 

1 + e. The estimates of Corollary 1.5 are uniform over all convex bodies C of all 
dimensions d. 
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The plan of the paper is as follows. In Section 2, we collect some facts needed 
for the proofs of Theorems 1.1 and 1.4. Namely, we review the classical result on 
the John decomposition of the identity operator and the minimum volume ellipsoid 
of a convex body, a recent result of Batson, Spielman and Srivastava [B+08] which 
allows one to obtain certain "sparsification" of the John decomposition, the stan- 
dard construction of tensor product from multilinear algebra which allows us to 
translate polynomial relations among vectors into linear identities among tensors 
and the classical construction of the Chebyshev polynomials which solve a relevant 
extremal problem. As it turns out, the vertices of the approximating polytopes P 
are picked up by certain algebraic conditions. 

We complete the proofs in Section 3. 

2. Preliminaries 

(2.1) Chebyshev polynomials. For a positive integer k let T k {t) be the Cheby- 
shev polynomial of degree k, see, for example, Section 2.1 of [BE95]. Thus for real 
t the polynomial (t) can be defined by 

Tk(t) = cos (/carccost) provided — 1 < t < 1 and 
T k (t) = \(t- \A 2 - l)^ + \( t + \A 2 - l)^ provided \t\ > 1. 
In particular, 

(2.1.1) |T fc (t)| < 1 provided \t\ < 1. 
Writing !&(£) in the standard monomial basis, we obtain 

m=0 v ' 

In particular, 

T x (t) = t, T 2 (t) = 2t 2 - 1, T 3 (t) = 4t 3 - 3t, T 4 (t) = 8t 4 - 8t 2 + 1. 

We note that Tk(—t) = Tk(t) if k is even and T/ c (— t) = —Tk(t) if k is odd. We also 
note that the polynomial Tfc(t) is strictly increasing for t > 1. 
In particular, 

(2.1.2) \T k (t)\ > ± ' K '- provided \t\ > r > 1. 

The polynomial Tj~(t) has the following extremal property relevant to us: for any 
to ^ [—1, 1] the maximum value of |p(to)|) where p is a polynomial of degp < k such 
that \p(t)\ < 1 for all t G [—1, 1], is attained for p = T k , see, for example, Section 
5.1 of [BE95]. 
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(2.2) Tensor power. Let V be Euclidean space with scalar product (•, •). For a 
positive integer k let 

v® k = V <g> • • • <g> V 



A; times 



be the /c-th tensor power of V. We consider V® k as Euclidean space endowed with 
scalar product (•, ■) such that 

k 

(xi <8> • • • <8> Xk, yi © ■ ■ ■ © Vk) = Y[{xi, Vi) 

i=l 

for all xi, . . . ,Xk,yi, . . . ,yfc G V. The space V® 2 is naturally identified with the 
space of all linear operators on V. 

The symmetric part Sym (V® k ) of V® k is the subspace spanned by the tensors 

x m = X ® • • • © X 
fc times 

for x £ V. The space Sym (V® fc ) is naturally identified with the space of all 
homogeneous polynomials p : V — > H. of degree fc. In particular, Sym (7® 2 ) can 
be identified with the space of quadratic forms on V and also with the space of all 
symmetric operators on V. We have 



dimSym(^) = ( dimy + fc - 1 ). 



Let us consider the direct sum 

W = R®V © 7® 2 © • • • © V m 

as Euclidean space with the standard scalar product, which we also denote by (•,•)• 
For a real univariate polynomial a(t) and a vector x 6 7, we denote by a® (a;) G 17 
the vector 

k 

(2.2.1) a®(x) = a © aix © a 2 x® 2 © . . . © a fc :r® fe , where a(t) = a m t m . 

m=0 

It is then easy to check that for any 6 7 and any polynomials a(t) and b(t), 
we have 

(a®(x), 6®(y)> = c((x,y», provided 
f2 2 2) ^ ^ k 

a(t) = Umt™, b(t) = Prnt™ and C(t) = ( a mf3 m ) t m . 

m=0 m=0 m=0 
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(2.3) The ellipsoid of the minimum volume. As is known, for any compact 
set C cM. d there is a unique ellipsoid of the minimum volume among all ellipsoids 
centered at the origin and containing C. If the minimum volume ellipsoid is the 
unit ball 

B = {xeR d : \\x\\ < 1}, 

where || • || is the Euclidean norm, the contact points Xi G CndB provide a certain 
decomposition of the identity operator i", called the John decomposition (recall that 
x eg) x for x G R d is interpreted as a d x d symmetric matrix) . We need the following 
result, see for example, [Ba97]. 

(2.3.1) Theorem. Let C C R. d be a compact set which spans ~§L d and let B C 
M. d be the unit ball. Suppose that C C B and that B has the smallest volume 
among all ellipsoids centered at the origin and containing C . Then there exist 
points xi, . . . , x n G C fl OB and non-negative real a\, . . . , a n such that 

n 

y^«i ® Xi) = /, 

1=1 

where I is the identity operator on R. d . Equivalently, 

n 

^2ai(xi,y} 2 = \\yf 
i=i 

for every y G M. d . 

(2.4) Sparsification. We need a recent result of Batson, Spielman and Srivastava 
on a certain "sparsification" of the conclusion of Theorem 2.3.1. Namely, we want 
to be able to choose the number n of points in Theorem 2.3.1 linear in the dimension 
d at the cost of a controlled corruption of the identity operator /. 

If A and B are dx d symmetric matrices we say that A ^ B if B — A is positive 
semidefinite. The following result is from [B+08] . 

(2.4.1) Theorem. Let 7 > 1 be a number and let x\, . . . ,x n be vectors in R d such 
that 

n 

Xi Xi = /, 

i=i 

or, equivalently, 

n 

^2(x z ,y) 2 = \\y\\ 2 
i=i 

for all y G M d . Then there is a subset J C {1, . . . , n} with \ J\ < "fd and (3j > for 
j G J such that 

jeJ 
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/ 7 + l + 2^7 \ 



or, equivalently, 



for all y G R d . 

3. Proofs 

We start with a lemma. 
(3.1) Lemma. Let C C M. d be a compact set. Then there is a subset X C C of 

\X\ < 4d 

points such that for any linear function £ : M. d — > R we have 

max\£(x)\ < max|£(:r)| < 3\fdmax.\£(x)\. 

xEX xEC xEX 

Proof. Without loss of generality we assume that C spans R d . Applying a linear 
transformation, if necessary, we may assume that C is contained in the unit ball B 
and that B is the minimum volume ellipsoid among all ellipsoids centered at the 
origin and containing C. By Theorem 2.3.1 there exist vectors xi, . . . x n e C n dB 
and numbers a±, . . . , a n > such that 

n 

y^«j (xi ®xi) = i. 

Applying Theorem 2.4.1 with 7 = 4 to vectors ^fcTiXi we conclude that for some 
J C {1, . . . , n} and /3j > for j e J we have 

(3.1.1) 7 =< >^n 7 ,•/) r< 9/ 
and |J| < 4d. We let 

X = { Xj : jeJ}. 

In particular, Xj G C and = 1 for all j G J. Comparing the traces of the 

operators in (3.1.1), we get 

(3.1.2) d < J2 a ifo - 9d - 

jeJ 



A linear function £ : R d — > R can be written as £(x) = (y, x) for some y G R d . It 
follows by (3.1.1) that 

J2( a jPj)(yi x j) 2 > \\y\\ 2 

jeJ 

and then by (3.1.2) it follows that 

\(y i x j)\ > 3^'^" for some j E J. 



Since C C B, we have 



max | (y, x) \ < 



and the proof follows. □ 

We now prove Theorem 1.1. 

(3.2) Proof of Theorem 1.1. Let us denote V = R d and let us consider the 
space 

W = R © V © V® 2 © • • • © V m , 
see Section 2.2. Let us define a continuous map <j> : V — > W by 

4>(x) = 1 © x © x® 2 © • • • © £® fc for x E V. 

We consider the compact set 

C = {(f>(x) : xeC}, C c W. 

We note that C lies in the subspace 

R © V © Sym (V® 2 ) © • • • © Sym (V® k ) . 

In particular, 

dunspan(c) < l + d+( 2 )+..-+( fc J = ( k )■ 
Applying Lemma 3.1 to C, we conclude that there is a set X C C such that 

such that for any linear function C : W — > R we have 

(3.2.1) max|£(0(x))| < max |£(0(x)) | < 3 max | 



We define P as the convex hull 

P = conv(XU-X). 

Clearly, P C C and P has at most 8( d ~£ fc ) vertices. To conclude that P approx- 
imates C with the desired accuracy, we compare the maxima of linear functions 
£ : R d — > R on C and on P. 
Suppose that 

£{x) = (y,x) for some y G V. 
Let us define a linear function C : W — > 1R by 

C(w) = (Tf(y), w) for all w e W, 

where T k is the Chebyshev polynomial of degree k, see Section 2.1 and (2.2.1). 
Then by (2.2.2), we have 

C(<f>{x)) =T k ({y,x)). 

Hence from (3.2.1) we obtain 

(3.2.2) mag\T k (£(x))\ < max \T k (£(x)) | < 3 ) max |T fc (^(x)) | . 

Suppose that £(x) < 1 for all x G P and hence \£{x)\ < 1 for all x G X. Then by 
(2.1.1) we have |T fc (£(o;))| < 1 for all x G X. If for some x G C we have £(x) > r 
then by (2.1.2) we have 

> (T _^^ (T + ^> ^ 

which contradicts (3.2.2). Therefore 

(3.2.3) max£(x) < max^(i) < rmax£(x) 

x6P x6C xEP 

for every linear function £ : R d — >• R, which proves that C C rP. □ 

(^3.3) Remark. One can sharpen the bounds somewhat by noticing that the polyno- 
mial T k is even for even k and odd for odd k. Consequently, the map (j) : V — > W 
can be replaced by 

<f> e (x) = 1 © x® 2 © • • • © x m ~ 2 © x m 

for even k and by 

4> (x) = x © x® 3 © • • • © x® fc " 2 © x m 

for odd /c. This allows us to replace ( d ~l~ k ) by D(d, k) defined by (1.1.1) throughout 
the statement of Theorem 1.1. 
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(3.4) Proof of Theorem 1.4. As in Section 3.2, we construct the space W, the 
map 0, the set (7 and the subset IcCso that (3.2.1) holds. We then define P as 
the convex hull 

P = conv(XU (-l/fj)X). 

Clearly, P C C and P has at most 8( d "£ fc ) vertices. To conclude that P approx- 
imates C with the desired accuracy, we compare the maxima of linear functions 
£ : R d — > R on C and on P. 

Let T k be the Chebyshev polynomial of degree k. We define a polynomial S k by 

S k (t) = T k (^-t+ fJ - 1 



// + 1 //+ 1 
Hence deg5 , / c (t) = fc. Moreover, 

(3.4.1) |£ fc (t)| < 1 provided - \i < t < 1 
and 

, (A-VA^I) fe +(A + yA^l) /c 

(3.4.2) |S fc (*)| > ^ '-— -i ^- provided * > r. 

Given a linear function £ : R d — > R, 

£(x) = (y, x) for some y G V, 
we define a linear function £ : W — > R by 

C(w) = (Sf(y), w) for all w G W. 

Then 

C(<t>(x)) = S k ((y,x)). 

Hence from (3.2.1) we obtain 

(3.4.3) max\S k (£(x))\ < max\S k (£(x))\ < Si J max |S fc (£(x)) | . 

Suppose that £(x) < 1 for all x £ P. Then, necessarily, 1 > £(x) > — ji for all x G X 
and hence by (3.4.1) we have \S k (£(x)) | < 1 for all x G X. If for some x G C we 
have £(x) > r then by (3.4.2) 

which contradicts (3.4.3). Hence (3.2.3) holds for every linear function £ : R d — > R 
and, therefore, C C rP. □ 
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(3.5) Proof of Corollary 1.2. Let us choose r = 1 + e in Theorem 1.1. We use 
the standard estimate 



Let us choose 

(3.5.2) k 
where /3 > is a constant. Then 



(3d 1 
-^ln- 



(3.5.3) ^(V^TI < ^(l + o(D), 



k y \ d J I - 4/3 

where "o(l)" stands for a term which converges to uniformly on d as e — > 0. 
On the other hand, 



In (t + Vt 2 - l) = v 7 ^ (1 + o(l)) , 



where "o(l)" stands for a term which converges to as e — )• 0. Then, as long as 

p> 1 



4V2 



the condition of Theorem 1.1 is satisfied for all sufficiently small < e < €q((3). 
The proof now follows by (3.5.1). □ 

(3.6) Proof of Corollary 1.5. To prove Part (1), we observe that A = 1 + 5 and 
hence 

(3.6.1) In (A + VA 2 - l) = V25 (1 + o(l)) , 

where "o(l)" stands for a term which converges to as 5 — > 0. In Theorem 1.4, 
let us choose k defined by (3.5.2) with e replaced by 5. Comparing (3.5.3) with e 
replaced by 5 and (3.6.1) we conclude the proof as in Section 3.5. 

To prove Part (2), in Theorem 1.4 we choose r = 1 + e and k defined by (3.5.2). 
Then 

/ \ 1/2 

(3.6.2) In (A + a/A 2 - l) = 2 {-^ J (l + o(l)), 

where "o(l)" stands for a term which converges to uniformly on \x > 1 as e — >■ 0. 
Comparing (3.6.2) and (3.5.3), we conclude that the condition of Theorem 1.4 is 
satisfied for all sufficiently small < e < e (/3) as l° n S as 



The proof now follows by (3.5.1). □ 
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(3.7) Proof of Corollary 1.3. Let us choose r = ^y/d/k in Theorem 1.1, where 
7 > is a constant. Using Stirling's formula, we conclude that for each k 



lim —6 lf , 

d >OC y/d \ 0, 



yf(i + «(D). 



where "o(l)" stands for a term which converges to as k grows. 
On the other hand, for each k 

r + y/r 2 - 1 27 
hm — = — . 

d >0O y/d y/k 

The proof now follows by Theorem 1.1. □ 
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